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£\\ Abstract. 

Photometric calibration is currently the dominant source of systematic uncertainty 
in exploiting type la supernovae to determine the nature of the dark energy. We review 
J> our ongoing program to address this calibration challenge by performing measurements 

of both the instrumental response function and the optical transmission function of the 
pH atmosphere. A key aspect of this approach is to complement standard star observa- 

O | tions by using NIST-calibrated photodiodes as a metrology foundation for optical flux 

measurements. We present our first attempt to assess photometric consistency between 
synthetic photometry and observations, by comparing predictions based on a NIST- 
diode-based determination of the PanSTARRS-1 instrumental response and empirical 
atmospheric transmission measurements, with fluxes we obtained from observing spec- 
i i trophotometric standards. 

> 

in 

^ 1. Introduction and Motivation 

The challenge of photometry is to extract knowledge of the location and flux distribu- 
tion of astronomical sources, based on measurements of the 2 dimensional distribution 

pvq of detected photons in a focal plane. Each pixel i in the detector array sees a signal S i 

y—\ given by 



)H sources j 



(A)Ri(A)T(A)G(A)Ai dA, (1) 



where the sum is taken over all sources (including the sky) that contribute to the flux in 
the pixel, $>j(A) is the photon spectrum for source j, Ri(A) is the throughput of the pixel, 
including the transmission of the optics and the pixel's quantum efficiency, T(A) is the 
optical transmission of the atmosphere, G(A) accounts for non-atmospheric extinction 
processes along the line of sight to the source (necessary to convert from top-of-the- 
atmosphere fluxes to the SED of the source), and A,- is the effective aperture of the 
system for pixel i, essentially the wavelength-independent part of the instrumental re- 
sponse. 

Photometric calibration uncertainties currently limit our ability to use type la su- 
pernovae to determine the nature of the dark energy ( [Sullivan et al.|20lT i. The essential 
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tasks are 1) to establish a clear relationship between the zeropoints in the various pass- 
bands used (to avoid contamination of the Hubble diagram with systematic zeropoint 
errors), and 2) to map out the shapes of the passbands (to enable precise K corrections). 
The absolute overall zeropoint of the system is degenerate with the intrinsic brightness 
of SN la's. Precise colors are what matter. However we do note that in splicing together 
a Hubble diagram from multiple surveys, it is important to understand the overall zero- 
points of the relevant instrumental photometric systems. 

There are a number of alternative approaches to establishing well- understood col- 
ors in a multiband photometric system, comprising some combination of the following: 

1. Using terrestrial sources to establish a spectral scale for celestial sources. This is 
the approach taken by Hayes & Latham|(|1975[) and Hayes et al. ( 1975 1, that in 



effect uses Vega as a transfer standard. In this instance the radiometric properties 
of the terrestrial source is used as the fundamental system calibration. 

Exploiting the theoretical understanding of stars to predict the SEDs of individ- 
ual stars. In practice white dwarfs are the favored class of object. Observations 
of appropriate sources are then used to calibrate instruments, and the theoreti- 
cal spectrum is the basis for the relative system calibration. This approach has 
been used for the HST calibrations described by |Colina & Bohlin ( |1994[ );[Bohlin 



(1996]) and presented at this meeting by Hunt. This approach is also being pur- 



sued by the SkyMapper project (as described by Bessell at this meeting). The 
status of theoretical spectral modeling is described by Rauch in his contribution 
to this meeting. 

Statistical properties of stars can also be used to establish consistency across the 
sky, and to correct for various types of chromatic attenuation. This approach was 
pursued by High et al. (2009 1, and w as used by|Schlafly & Finkbeiner| ( 2011[ ) to 



show that the commonly-used SFD (Schleg eTet al.||1998 ) extinction map (G(A) 



of equation 1) requires renormalization of 0.86 in E{B - V) compared to SFD. 



4. An alternative approach, described in Stubbs & Tonry| (|2006[) and in the contri- 



bution by Cramer to this meeting, attempts to bring modern metrology methods 
to bear on the calibration challenge. This approach uses NIST-calibrated detec- 
tors ( Ahtee et al.|2007 1 as the fundamental reference for establishing the system's 



sensitivity function. Initial results from this approach were obtained on the CTIO 
4 meter (Stubbs e t al.||2007b| ), and more recently on the PanSTARRS-1 system 



( Stubbs et al.||2010 l. The experimental challenge in this approach is to obtain a 
reliable measurement of the instrumental response function that is traceable to SI 
standards. 

We stress that these approaches can all be pursued in parallel, and the relative 
consistency we achieve can be used to understand, quantify, and overcome sources of 
systematic error. 

2. Throughput Measurements of the PanSTARRS-1 Survey System 

2.1. Full-aperture system response function measurements 

We have used a photodiode to monitor the flux emanating from a back-illuminated flat- 
field screen in the dome of the PanSTARRS-1 telescope to map out the full-aperture 
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relative system response function. Details of this measurement are given in Stubb s et al. 
( 20 1 ) and Tonry et al. ( 2012[ ) . The results are shown in Figure [T] We project light from 
a tunable laser onto the flat-field screen in the dome. We measure the flux emanating 
from the screen, incident on the telescope pupil, with a calibrated photodiode. We then 
compare the flux detected by the instrument to the incident flux, as measured by the 
photodiode. Performing this measurement at a succession of wavelengths allows us to 
determine system throughput as a function of wavelength, using the calibrated photo- 
diode as the fundamental reference. A philosophically similar full-aperture calibration 
system for the Dark Energy Survey is described in these proceedings by Marshall. 
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Figure 1 . Shown on the left is the measured relative system response function for 
PS-1, all normalized to a NIST photodiode. The panel on the right shows the radial 
position-dependence of the filter transmission function. Each curve corresponds to a 
different radius on the focal plane. 



2.2. Collimated-Light Determination of Ghosting in the PanSTARRS-1 System 

During the processing of the throughput data described above, we became increasingly 
concerned about stray light, multipath effects, and ghosting in the optical train. The 
basic problem is that by imaging a uniform surface brightness screen we can't distin- 
guish the focussing light paths through the system (which is how we measure celestial 
sources) from other light paths. This is the underlying reason why dome flats, twilight 
flats and sky flats typically mutually disagree, and require an "illumination correction". 
Especially at wavelengths where the filter transmission is low, we find a substantial 
amount of light scatters from the focal plane, up to the filter, and back down to the 
CCDs. This is a source of systematic error. The contribution by Regnault to this meet- 
ing describes the merits of ray tracing to model these ghosts, and the LSST calibration 
team is undertaking a similar exercise (see paper here by Lynne Jones). 

In order to quantify the amount of ghosting in the system, we set up a collimated 
telescope that was fed by an optical fiber to send a beam of collimated light onto the 
primary. This controls the phase space distribution of rays entering the telescope much 
better than the light emanating from the flat field screen. Moreover, we can adjust the 
focus of the collimating telescope to control the size of the image on the PS-1 focal 
plane. Figure [2] shows an example image obtained with this configuration. 

We used the ghost-light fraction as a function of wavelength to correct the flat- 



screen generated throughput data. This is described in more detail in |Tonry et al. ( 2012 1. 



4 



Stubbs, Tonry 




Figure 2. This is a false-color image of the light distribution across the entire 
PanSTARRS-1 focal plane when illuminated by a collimated beam. The lack of 
axisymmetry in the beam spot is because the calibration diode obscured the output 
pupil of the collimation telescope. We used this determination of the scattered light 
vs. wavelength to obtain a more reliable measurement of relative system throughput. 



3. Atmospheric Transmission 



Our team's review paper on the atmosphere (Stubbs et al. 2007a) describes the vari- 
able components of atmospheric transmission that require attention in order to achieve 
improved photometric performance. Water vapor, ozone, clouds and aerosols are the 
primary concerns. Rayleigh scattering from 02,-^2 and other well-mixed gases is es- 
sentially deterministic, are the molecular absorption line strengths from Oj- Although 
a number of groups are conducting or are planning to undertake explicit atmospheric 
transmission monitoring (see contributions to this conference by McGraw, the Texas 
A&M team, the ESO program, and Blake), it remains the case that our community 
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does not yet know the angular and temporal correlation functions of these variable at- 
mospheric attenuation processes. Max Fagin and Justin Albert's presentations to this 
conference describe our group's program to fly a set of laser diodes to characterize 
aerosol extinction, so we will focus here on our measurements of water vapor. 

3.1. Water Vapor Measurements at the Pan-STARRSl site, using Polaris 

We have constructed an objective grating slitless imaging spectrograph (Shivvers et 
al. in prep). The configuration is very similar to the spectrophotometric system de- 
scribed by McGraw. We decided to point at the North Celestial Pole, so as to obtain 
data at fixed airmass over an extended period of time. The camera in this system uses 
a Pixis 1024BR deep-depletion detector and so it has low fringing. The dispersion 
and focal length are such that we can obtain a spectrum spanning 300 to 1000 nm, 
and the 50-50 objective transmission grating greatly suppresses second order contam- 
ination light. Figure [3] shows the comparison between one of our 1-d spectra and the 
MODTRAN model that has been tuned to obtain an equivalent width of water vapor ab- 
sorption to match the observations. The instrument was installed on Haleakala in July 
2011 and has been acquiring data intermittently ever since. Figure [4] shows the tempo- 
ral evolution, over a one month period, of the equivalent widths (EW, in nm) of the O2 
and water vapor features. We obtain an excellent match to the MODTRAN prediction 
of the EW for O2, and there is no evidence of this feature changing over the duration 
of the observations. Both of these facts indicate that we are obtaining meaningful and 
reliable measurements of absorption features. The water vapor attenuation does show 
considerable variation over this period, and we assess our fractional accuracy in PWV 
determination to be about 10%. This is adequate to obtain an overall precision of < 1% 
in the determination of transmission in the Pan-STARRSl passbands. 

3.2. MODTRAN, tweaked 

For comparison with the standard star observations, we used the MODTRAN "Generic 
Tropical" model atmosphere, with the "Desert Extinction (Spring-Summer)" aerosol 
choice. No attenuation from clouds was included. The PWV at the base of the at- 
mosphere was set to 0.65 cm, since this produced a good match between our Polaris 
observations on the night the standard stars were observed, and the MODTRAN predic- 
tion for the line of sight from the altitude of the summit of Haleakala, looking towards 
the celestial pole. As described below we ended up making an adjustment to the aerosol 
component, based on the airmass dependence we observed. We note that the other, de- 
terministic, aspects of atmospheric transmission were not adjusted. We used the default 
ozone column. This amounted to a MODTRAN transmission model with three empiri- 
cally adjusted parameters: PWV, aerosol optical depth, aerosol Angstrom exponent. 



4. The Total PanSTARRS-1 System Throughput: Instrument and Atmosphere 

Taking the ghosting-corrected relative system throughput function in conjunction with 
the MODTRAN atmosphere (with parameters adjusted based on our observed PWV), 
we obtained the PanSTARRS system throughput functions shown in Figure [5] We have 
adopted a convention of 1.2 airmasses as the definition of the PS throughput function 
( |Tonry et al.|2012] |. 
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Figure 3. This Figure compares the atmospheric transmission measured with the 
objective grating dispersed imager with the model atmosphere produced by MOD- 
TRAN, once we tuned the precipitable water vapor to match the equivalent widths 
seen in the data. Note the Polaris spectrum also has features from the stellar atmo- 
sphere, such as the Ca triplet at 850 nm. 



With the PS-1 system throughput function in hand, we have used measurements 
of HST spectrophotometric standards to compare our NIST-based flux calibration to 
the fluxes observed from these standards. Details of this are presented in |Tonry et al. 



(2012). An important aspect of this comparison is to determine the atmospheric trans- 
mission. 

For the Pan-STARRS 1 bandpasses we integrated a set of power law SEDs against 
each of these model atmospheres for each bandpass and created an interpolation func- 
tion for the extinction as a function of four variables: z for airmass (sec£ where £ is 
the zenith distance), h for precipitable water vapor (PWV) (typically 0.65 cm at sea 
level), a for "aerosol exponent" (nominally 1; we modify the Modtran aerosol compo- 
nent by applying this power to the transmission, thereby mostly affecting the aerosol 
amplitude), and p for SED power law (we calculate the extinction for pure power law 
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Figure 4. This Figure shows the evolution of the measured equivalent widths of 
O2 and H2O over the duration of one month. Both the accuracy and the precision of 
the oxygen line indicate that we are making a reliable measurement of atmospheric 
attenuation. 



SEDs, where p - +2 for f v ~ v corresponds to an O star with (r—i) - -0.43, p = 
for f v ~ const corresponds to an F star with (r—i) = 0.00, and p = -2 for f v ~ v~ 2 
corresponds to an K5 star with (r—i) = +0.42. The extinction dm in magnitudes is 
parameterized the extinction as 

\ndm - InC + Zlnz + A In a + Pp + \nh(H + H\ \nz + H 2 \nh) (2) 

The coefficients for each of the Pan-STARRSl filters are given in Table [T] 

5. Standard Star Observations 

MJD 55744 (UT 02 July 2011) was a photometric night during which we observed a 
substantial number of spectrophotometric standard stars from the STIS Calspec (Bohlin 
et al. 12001) tabulation: 1740346, KF01T5, KF06T2, KF08T3, LDS749B, P177D, and 
WD 1657- 343. These were observed throughout the night at airmasses between 1 and 
2.2 in all six filters and also with no filter in the beam. Each observation was repeated, 
and exposure times were chosen to stay well clear of any non-linearities but still permit 
good accuracy. Observations in ypi of the fainter white dwarfs were curtailed at 100 
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Figure 5. Combining the ghost-corrected instrumental response function with the 
PWV-informed MODTRAN atmosphere at 1 .2 airmasses produces the system sen- 
sitivity functions shown here. 



sec duration, so their uncertainties are relatively large. In addition, Medium Deep Field 
09 (which overlaps SDSS Stripe 82) was observed a dozen times in each of gp\ rp\ ip\ 
zpi and ypi offering the opportunity to tie the spectrophometric data to a well-observed 
Pan-STARRS 1 field. All stars were placed on OTA 34 and cell 33, so their integration 
was on the same silicon and used the same amplifier for read out (gain measured to 
be 0.97 e~/ADU). The observations were bias subtracted and flatfielded as part of the 
normal IPP processing, and the IPP fluxes (instrumental magnitudes) were then avail- 
able for comparison with tabulated SEDs. The IPP performs an aperture correction and 
reports fluxes within a radius of 25 pixels (13 arcsec diameter). Observations of Polaris 
on MJD 55744 with the spectroscopic sky probe had a PWV indistinguishable from the 
long term mean of 0.65 cm. 



6. Putting it all Together: Closing the Photometric Loop and Assessing Consis- 
tency 



As described in more detail in Tonry et al. (2012 1, we compared the observed gp\ rp\ 



ip\ zp\ and yp\ fluxes for the HST spectrophotometric standards with the predictions 
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Table 1. PanSTARRS extinction coefficients. The columns contain coefficients 
described above for each of the PanSTARRS bandpasses that interpolate the Mod- 
tran extinction calculations. The final column is the percentage scatter of these fits 
relative to the calculated values. Note that the saturation of the water lines means 
that the extinction is not proportional to sec f (Z + 1), particularly for ypi. This is 
consistent with the modeling results from MODTRAN. 
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obtained from using equation (1), using our in-dome determination of the instrumental 
response, the observationally adjusted MODTRAN model for the atmosphere, and the 
HST Calspec data for the source SEDs. We find that we need to apply a gently varying 
"tweak" to the system response function, with an rms of around 3-4%, in order to 
obtain consistency between the synthetic photometry using CalSpec SEDs and the on- 
sky observations. 

7. Conclusions 

The results presented here comprise, to our knowledge, the first instance of a NIST- 
calibrated telescope response function and a MODTRAN atmospheric model gener- 
ating synthetic photometry that is then compared with on-sky measurements of spec- 
trophotometric standards. We obtain agreement at the 5% level, except for the ypi band 
where the discrepancy is 10%. We consider this to be encouraging, since it's our first 
attempt to "close" this photometric loop. Our objectives for future work include: 

1. Obtaining more precise on-site determination of atmospheric aerosols. 

2. Improving the in-dome calibration technique, and our corrections for optical mul- 
tipath effects. 

3. Identifying the origin of the "tweak" we need to apply to obtain overall con- 
sistency between the NIST-based calibration and the expectations from stellar 
SEDs. 
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Discussion 

Q: (S. Deustua) What about flatfields? How reliable is your flat/illumination system in 
that a source placed on any pixel gives the same result? You have a very wide (3 deg) 
field. 

Q: (N. Regnault) It is difficult to apply a photometric calibration obtained from a flat- 
field-like illumination to point source photometry. Do you have specific plans to bridge 
this gap? 

A: We agree that it's difficult to distinguish the focusing light paths from the stray and 
scattered light that lands on the focal plane, when the illumination comes from a dif- 
fuse flat field screen. As described above, we have used a collimated beam to attempt 
to disentangle these illumination paths, and to determine the response function at one 
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location on the CCD array. The task of tying together the system sensitivity across the 
array is a distinct but related problem. A number of groups are working to develop var- 
ious calibration devices that control the phase space distribution of the rays (in angle 
and position) that are sent into the telescope pupil from a monitored calibration source. 
We think that making an "illumination correction" is an essential ingredient in reduc- 
ing and assessing sources of systematic error. This is best obtained on the sky, either 
by rastering stable sources across the focal plane, or through the "ubercal" procedure 
that uses the multiple observations per field from a wide-field multi-epoch survey. As a 
community we are still developing an overall methodology for obtaining the best possi- 
ble accuracy and precision, and this meeting has been very useful in helping move this 
forward. 



